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Abstract. There is a remarkable and canonical problem in 3D geometry and 
topology: To understand existing models of 3D fluid motion or to create new 
ones that may be useful. We discuss from an algebraic viewpoint the PDE 
called Euler's equation for incompressible frictionless fluid motion. In part I 
we define a "finite dimensional 3D fluid algebra," write its Euler equation and 
derive properties related to energy, helicity, transport of vorticity and linking 
that characterize this equation. This is directly motivated by the infinite di- 
mensional fluid algebra associated to a closed riemannian three manifold whose 
Euler equation as denned above is the Euler PDE of fluid motion. The clas- 
sical infinite dimensional fluid algebra satisfies an additional identity related 
to the Jacobi identity for the lie bracket of vector fields. In part II we discuss 
informally how this Jacobi identity can be reestablished in finite dimensional 
approximations as a Lie infinity algebra. The main point of a developed version 
of this theory would be a coherence between various levels of approximation. 
It is hoped that a better understanding of the meaning of the Euler equation 
in terms of such infinity structures would yield algorithms of computation that 
work well for conceptual reasons. 



1. Algebra and Topology of Ideal Fluids 

A finite dimensional 3D-fluid algebra is a finite dimensional vector space V 
provided with three structures: 

1) an alternating trilinear form { , , } on V, called the triple intersection form. 

2) a symmetric nondegenerate bilinear form ( , } on V , called the vorticity 
linking form. 

3) a positive definite inner product ( , ) on V", called the metric. 

If M is a 3D closed oriented Riemannian manifold there is a classical example 
of a fluid algebra which is infinite dimensional and which is constructed inside the 
differential forms: 

• V consists of the coexact one forms, the image of two forms by the operator 
*d* where * is the Hodge star operator of the metric. Under the correspon- 
dence between one forms and vector fields given by the metric, elements in 
V correspond to volume preserving vector fields which have flux or net flow 
zero across any closed hypersurface. 

• The vorticity linking form on V is defined by setting 
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This is equal to J M (da) A b, depends only on da and db, and may be con- 
strued as a linking number of the two one-dimensional transversally mea- 
sured foliations defined by the kernels of da and db and transversally mea- 
sured by the two forms da and db (see Arnold-KhesinjAK and Sullivan|S]). 
Here da and db determine the vorticities of the vector fields corresponding 
to a and to b. da and db may be approximated, in the sense of integrating 
against a smooth test form, by weighted sums of closed curves which bound 
weighted sums of surfaces approximating a and b. Recall that linking num- 
bers are defined by intersecting one set of curves with surfaces bounding 
a second set of curves. Thus the integral approximately computes the to- 
tal linking of the weighted family of curves approximating da with those 
approximating db. 

One can see from the Hodge decomposition that the vorticity linking 
form is nondegenerate on coexact one-forms. 
• The triple intersection form {a, b, c} is the integral of a wedge b wedge c 
over M, which may be construed as a triple intersection of the surfaces with 
boundary approximating a, b, and c. 

The Euler evolution in V describing ideal incompressible frictionless fluid mo- 
tion in M is an ODE whose solutions may be described in words as follows: an 
isotopy from the identity is a fluid motion iff the path in V corresponding to the 
instantaneous velocity of the motion satisfies: the vorticity at time t is the two form 
transported by the motion from the vorticity two form at time zero. 

We can write out this Euler ODE just using the elements of a fluid algebra. 

Namely, the right hand side of the evolution ODE for X(t) is described by its 
inner product with any vector Z as follows: (^jr, Zj = {X, DX, Z} where D is the 
operator on V defined by (DX,Y) — (X,Y). D is called the curl operator on the 
elements of a fluid algebra. Note that (DX, Y) = (X, DY) since ( , ) is symmetric 
by assumption. 

Proposition 1.1. The "energy" — (X, X) and the "helicity" — (X,DX) are each 
preserved by the evolution in a fluid algebra. 

Proof. 1) ^(energy) = 2(^-,X) = 2{X,DX,X} = by the alternating 
property of { , , }. 

2) -| (helicity) = 2 DX) = 2{X, DX, DX} = again by the alternating 
property of { , , }. 

□ 

Corollary 1.2. For a finite dimensional fluid algebra, there is no finite time blowup 
and the flow stays on the intersection of the energy sphere with the helicity level 
set. 

Notice in this classical three manifold case the alternating trilinear form on all 
one-forms and the bilinear form on all one-forms are purely topological and don't 
depend on the metric. The metric is needed to define the subspace of coexact forms 
inside all one-forms to which one restricts the trilinear form, the bilinear form and 
the inner product to obtain the infinite dimensional fluid algebra. 

Also notice that even though the Proposition follows for this classical Euler case 
by the same short proof, the Corollary of course does not follow because of the non 
compactness of the infinite dimensional sphere. The truth or falsity of infinite time 
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existence of smooth solutions of the Euler evolution for arbitrary smooth initial 
conditions is a celebrated open problem. 
Here is another family of examples: 

Example 1.3. Given a finite dimensional Lie algebra [, ] with a nondegenerate 
invariant symmetric pairing ( , ) define a fluid algebra by: 

• for the linking form take ( , ) . 

• for the intersection form take {X, Y, Z} =([X, Y], Z) and 

• for the metric take any positive definite inner product. 

Euler flows for these examples were described differently by Arnold and Khesin 
[AKj . The interested reader may verify that their evolution equation (see |AK[ 
Theorem 1.4.9]) is related to ours by a linear coordinate change. 

Notice: 

1) For a general fluid algebra one may reverse the above formula to define a 
bracket by the formula ([X, Y],Z)= {X, Y, Z}. If the Jacobi identity were 
satisfied for this bracket the fluid algebra would arise from the Lie algebra 
example. 

2) In the infinite dimensional example related to 3D manifolds and the Euler 
ODE this Jacobi identity is satisfied and corresponds to the Lie algebra of 
volume preserving flux zero vector fields. 

3) In order to have models which satisfy the invariance of energy, helicity, 
and more generally circulation (see below) and the transport property of 
vorticity (see below) this Jacobi relation is not required. The significance 
of the Jacobi relation needs exploring (see Part II). One interpretation 
of Jacobi is the following: if we extend the transport (defined below) to 
higher tensors by the derivation property then Jacobi is equivalent to the 
statement that transport fixes the tensor defining transport itself. 

Definition 1.4. (Velocity and Vorticity) We think of elements in V as analogous to 
velocity and vorticity fields of the ideal fluid. The vorticity of a velocity is obtained 
by applying the curl operator D to the velocity. 

Proposition 1.5. If X(t) satisfies the Euler equation of the fluid algebra 




= {X,DX,Y}, 



then the vorticity DX(t) = Y(t) satisfies 

(^)= { .™ 

where D' is the inverse of D. 
Proof. j,DX(t) = D(±X(t)). So 

= {X, DX, DZ} = {X, V, DZ} = {D'Y, Y, DZ}. 

□ 
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Definition 1.6. (Transport) Define the infinitesimal transport of Z by X, T(X, Z), 
by the condition that its inner product with arbitrary W is given by (T(X, Z), W) = 
{X,Z,DW}. 

Note T(X, Z) = -T(Z, X), but Jacobi is not assured for T(X, Y). The transport 
is meant to model for volume preserving vector fields the Lie derivative action, which 
by the way is canonically equivalent to the the Lie derivative action on closed two 
forms. 

Corollary 1.7. // X(t) satisfies the Euler equation of a fluid algebra then the 
vorticity Y(t) = curlX(t) = DX(t) satisfies dY Jf^ = T(X(t), Y(t)), namely the 
vorticity is transported from one time to another by the motion. 

Proof. Since Y = DX, by the proposition and the definition of transport the inner 
product of each with Z equals {X, DX, DZ}. □ 

Note: This definition of transport may be derived in the classical case from the 
expression of the Lie bracket of two volume preserving vector fields [V, V] in terms 
of the corresponding one-forms by the formula [V, V] = *d * (V A V) where *d* is 
the adjoint of d. 

Corollary 1.8. (Invariance of circulation) If the vorticity DZ(t) of a time depen- 
dent field Z(t) satisfies -^(DZ(t)) = T(X(t), DZ(t)), then the linking number of 
DZ(t) and DX(t) stays constant. 

Here, consistently with the situation discussed in the classical example, the link- 
ing number of DZ(t) and DX(t) refers to the value (X, Z) of the vorticity linking 
form on the fields X and Z . 

Proof. The linking number of DZ and DX is equal to (X, DZ). So 

l { x,DZ) = ^ t x(t),DZ(t)) + (x,±DZ { t)) 

= {X, DX, DZ} + (T(X, DZ),X) 
= {X, DX, DZ} + {X, DZ, DX} 
which is zero by the alternating property. □ 

Conclusion 1.9. The Euler ODE for velocity associated to a fluid algebra is the 
only ODE whose evolution keeps constant the linking of the transported vorticity of 
a general field with the evolving vorticity of the velocity solution of the ODE ( as 
expressed in Corollary \1.8\) . 

Proof. The calculation of the proof of Corollary 11.81 shows cancellation takes place 
iff the evolution of the velocity X(t) satisfies ( f t X(t), DZ) = {X, DX, DZ}. Since 
D is invertible this is equivalent to the definition of the Euler ODE. □ 

2. Algebraic Topology of Ideal Fluids 

In the classical example of ideal fluid motion on a metric three manifold M 
without boundary, the fluid algebra is embedded in the differential algebra of all 
differential forms. So we are adding to the fluid algebra non constant functions, 
non constant volume forms, more one-forms and more two forms. New operators 
and operations appear: exterior derivative, the wedge product, the integral of three 
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forms over the manifold, the Hodge star operator and the Hodge decomposition with 
the associated projections: all forms = exact forms + coexact forms + harmonic 
forms, (exact means image d and coexact means image of *d*). The fluid algebra 
in these terms is made from the subspace of coexact forms in degree one, the curl 
operator obtained by doing d then star, the alternating form and metric as described 
above using wedge, star and the integral. Alas, all of this is infinite dimensional so 
compactness and the easy proof above of long time existence disappears. 

However there are many finite dimensional models of the algebraic topology 
residing in the above structure. Two types are: 

I (grid type). Divide the manifold into cells and average the forms over 
oriented faces of the correct dimension. By Stokes' theorem one gets a map 
of chain complexes 7:(all smooth forms with d) — »■ (cellular cochains with 
d) which induces an isomorphism on cohomology. There are maps of chain 
complexes in the opposite direction using the heat flow and dual cells which 
are inverse to / up to chain homotopy and invariant under star (from class 
lectures and an unpublished manuscript). 
II (eigenvalue type). Consider the eigendecomposition of all p forms for the 
laplacian operator (— l) p ~ 1 (*d*d— d*d*) and project onto finite dimensional 
pieces by putting bounds on the eigenvalues. These are invariant under d 
and star as well and induce isomorphisms on cohomology. For the flat torus 
these are used in numerical calculation for fluids. 

A third type which is due to Whitney, is very elegant and useful for algebraic 
topology (see |S3j ) but it is not invariant under star for a fundamental reason which 
to my knowledge has never been corrected. So we do not consider Whitney forms 
here. 

The first two finite dimensional approximations have everything they need for 
defining the fluid algebra approximating the classical fluid algebra for the classi- 
cal Euler evolution EXCEPT the wedge product. This lack can be corrected by 
forming the wedge product on the finite dimensional image inside forms and then 
integrating or projecting back to the finite dimensional model. This is how we mul- 
tiply numbers on a computer using finite parts of the decimal expansion and how 
numerical computations for fluids are performed. In each case the multiplication 
obtained is appropriately commutative but NOT associative. 

Now nontrivial ideas of algebraic topology enter. First recall that a chain map- 
ping of chain complexes of vector spaces (or free Z-modules) inducing a homology 
isomorphism has an inverse up to chain homotopy. Such inverses and the chain ho- 
motopies can be used to transport algebraic structures (up to homotopy) between 
chain complexes that are very different as vector spaces, for example the infinite 
dimensional deRham complex and the finite dimensional models. The associator 
of this finite dimensional wedge product ((a A b) A c) - (a A (b A c)) is a three to 
one operation which is a mapping of chain complexes which commutes with the 
the natural differentials. Moreover because of the homology isomorphism above 
this associator is the commutator with d of a correcting three to one operation of 
degree one less which we denote Ai. Now we continue using ideas of Stasheff's 
Princeton Thesis (1959). Consider the five ways to associate four entities. These 
arrange naturally at the vertices of a pentagon. We sum the corrections Ai on each 
edge combined with A to build a four to one operation of degree — 1 which again is 
a cycle in that it commutes with d extended to multivariables by the Leibniz rule. 
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By the homology isomorphism above there is a degree —2 four to one operation 
A2 which fills in the cycle corresponding to the boundary of the pentagon, i.e. one 
whose commutator with d is the four to one operation of degree —1 obtained by 
going around the pentagon. This process continues indefinitely producing n to 1 
correction operations of degree — n whose commutator with d is an appropriate 
more and more complex looking formula in the inductively constructed corrections. 

There are two versions of this: the original Stasheff one is very beautiful. The 
complex looking formulae are nothing more than the combinatorics of the moduli 
spaces of complete hyperbolic surfaces with geodesic boundary boundedly related 
to the two disk with three or more punctures on the boundary, then naturally 
compactified by geometric limits. This is the model also controlling genus zero part 
of open string theories. The Stasheff polyhedra can also be described as a moduli 
space of planar rooted trees. 

The second way remembers the commutativity and seems more appropriate here. 
It may be modeled on the moduli spaces of rooted trees (in space) whose leaves 
are labeled and whose interior edges are painted black or white. This model is a 
picture of what is called the bar cobar construction for algebras over an operad (see 
VallettejVj and WilsonfW]). 

Interestingly enough there is a specific procedure for computing these corrections 
by placing the chain homotopies mentioned above on the interior edges of the trees 
and using the wedge at the vertices. This procedure is identical to the tree part 
of the Feynman diagram algorithm in perturbative Chern-Simons quantum field 
theory where the "propagator" there is the chain homotopy here. 

After doing this work we obtain the derived or "infinity" version of the graded 
commutative wedge product compressed onto the finite dimensional approximating 
models. 

The same ideas may be applied to the bracket or transport discussed above 
which satisfies Jacobi in the infinite dimensional model. The Jacobi identity may 
be encoded in the the algebraic statement that on forms the wedge product and 
the adjoint of d = (— l) p * d* satisfy the following identity in words: the deviation 
of *d* from being a derivation of A is, as a two variable operator, itself a derivation 
in each variable. Adding the equation that *d * *d* = yields formally that the 
two variable operation satisfies Jacobi. This Lie bracket on forms becomes via 
the metric isomorphism the usual Lie bracket of vector fields extended by Leibniz 
to all multivector fields. (This is referred to as the Schouten-Nijenhuis bracket.) 
This relationship between *d* and A means we have a so called Batalin-Vilkovisky 
algebra or briefly a BV algebra. The BV formalism in perturbative quantum field 
theory due to Batalin-Vilkovisky is perhaps the most natural for mathematicians 
(see Costello jC]). The BV structure (*d*,A) may also be compressed to a new 
structure of BV algebra up to homotopy into our finite dimensional models. The 
bracket so obtained is called a Lie infinity structure. 

So now we have a commutative infinity structure related to d and a Lie infinity 
structure related to *d* whose leading term is the transport above. 

Before rushing off to make models that would be used for fluid simulation these 
ideas need I think to be completed in the following way: 

1) describe the complete structure used above d,star, A and the integral as 
an algebraic structure. Then develop the diagrammatic compression algo- 
rithms for this structure which yield the derived or infinity version of the 
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structure on the finite dimensional approximations. To my knowledge the 
current abstract homotopical algebra just falls short of this task. General 
multilinear operations with outputs can be treated as in |S2| . but pairings 
are not treated there or anywhere else to my knowledge (although Kevin 
Costello's work on renormalization and perturbative quantum field theory 
comes close [C]). 

2) understand the Euler evolution as a functorial construction on the derived 
or infinity version of this algebraic structure. The efforts of Arnold and 
Khesin and others are a beginning but to my knowledge this goal is not yet 
achieved (although describing the flow by an action principle may lead to 
such a functorial principle). 
If these two tasks are completed, we will have conceptually natural effective theories 
of fluid motion at every scale that fit together in an appropriate sense. There will 
be, by definition almost, natural algorithms for fluid computation based on the 
corrections that emerge from the compression of the algebraic structure into finite 
dimensions. These may play a role in proving long term existence of the classical 
ODE if that long term existence is true. If it is not true, this may also be revealed 
in these models which work at every scale. In either case we will have natural 
algorithms for computations of real fluids at every scale which of course do have 
long time existence and which are potentially observable at extremely small scales. 
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